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Abstract 

During the last two years, Europe has been facing a debt crisis, and Greece has been at its center. 
In response to the crisis, drastic actions have been taken, including the halving of Greek debt. Policy 
makers acted because interest rates for sovereign debt increased dramatically. High interest rates 
imply that default is likely due to economic conditions. High interest rates also increase the cost 
of borrowing and thus cause default to be likely. In equilibrium markets, economic conditions are 
used by the market participants to determine default risk and interest rates, and these statements 
are mutually compatible. If there is a departure from equilibrium, increasing interest rates may 
contribute to — rather than be caused by — default risk. Here we build a quantitative equilibrium 
model of sovereign default risk that, for the first time, is able to determine if markets are consistently 
set by economic conditions. We show that over a period of more than ten years from 2001 to 2012, 
the annually-averaged long-term interest rates of Greek debt are quantitatively related to the ratio 
of debt to GDP. The relationship shows that the market consistently expects default to occur if 
the Greek debt reaches twice the GDP. Our analysis does not preclude non-equilibrium increases 
in interest rates over shorter timeframes. We find evidence of such non-equilibrium fluctuations in 
a separate analysis. According to the equilibrium model, the date by which a half-default must 
occur is March 2013, almost one year after the actual debt write-down. Any acceleration of default 
by non-equilibrium fluctuations is significant for national and international interventions. The 
need for austerity or other measures and bailout costs would be reduced if market regulations were 
implemented to increase market stability to prevent the short term interest rate increases that 
make country borrowing more difficult. We similarly evaluate the timing of projected defaults 
without interventions for Portugal, Ireland, Spain and Italy to be March 2013, April 2014, May 
2014, and July 2016, respectively. The markets consistently assign a country specific debt to GDP 
ratio at which default is expected. All defaults are mitigated by planned interventions. 
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I. MOTIVATION 



Europe has been facing what the media has called "a ferocious debt crisis" since the 
beginning of 2010 [1]. Several European governments have accumulated what many consider 
to be unsustainable levels of government debt [2]. Greece has been at the center of the 
crisis, with the highest levels of public debt in the Eurozone and one of the biggest budget 
deficits [3]. Fears of default and its consequences [4] have led to special bailout funds [5], 
austerity programs in Greece [6], and recently agreements by banks to voluntarily dismiss 
half of Greece's debt [7]. Policy makers have taken action because of the dramatic increase in 
interest rates, which curtails the ability to borrow money from international capital markets. 
High interest rates imply that the market considers the risk of default to be high due to 
economic conditions, but they also increase the cost of borrowing and cause default to be 
more likely. Because interest rates both reflect expected default risk and cause default risk, 
it is not entirely clear which is the actual cause of default: economic conditions, i.e. an 
equilibrium market, or non-equilibrium market behavior. 

In equilibrium economics, market prices are based upon value and value is determined 
by price. This self-consistency in equilibrium hides the important subtle causal loops of 
the interaction of price with value that become manifest when deviations from equilibrium 
occur. Among the causes of deviations from equilibrium are herd behaviors that cause 
price changes whose timing is not related to inherent changes in value, but only to the 
intrinsic dynamics of the market itself. Such deviations from equilibrium may play a role 
in fluctuations characterized by "volatility" (though volatility can also include variations in 
equilibrium prices and value) but the impact of fluctuations away from equilibrium is poorly 
understood because they depend on those poorly characterized causal loops. 

In recent years, major market indices have varied widely, amid concerns about macroeco- 
nomic conditions that affect broadly the values of goods and services. In the case of markets 
whose values depends on macroeconomic conditions, the equilibrium market price reflects 
economic news. When market price is itself a causal factor in economic conditions, which 
it often is, the feedback to economic conditions can lead to inherent dynamics that are not 
found in the theory of equilibrium prices. Because markets are considered to represent eco- 
nomic conditions, the interpretation of market changes is used in policy setting for economic 
interventions. When non-equilibrium fluctuations of the market are interpreted as indica- 
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tors of economic conditions for policy interventions, this misinformation can lead to policy 
actions that are not justified by economic conditions but only by the false signals given 
by market fluctuations. However, in some markets, prices directly impact on the economic 
system. This is the case in bond markets, where interest rates that are set by the market 
directly impact on the borrowing costs of countries. Whether the markets are or are not in 
equilibrium is therefore key to both economic conditions and correct policy response. Disen- 
tangling the perception of market fluctuations from economic conditions, and the possibility 
of market fluctuations determining rather than reflecting economic conditions is important 
to both economic theory and policy. While equilibrium prices are important components of 
the proper and constructive function of economic markets, the impact of non-equilibrium 
fluctuations may be destructive, and therefore should be the subject of policy attention, 
which can by careful choice limit both the size and impact of non-equilibrium volatility. 

Here we build a simple equilibrium model of sovereign default risk that, for the first 
time, is able to directly quantify the consistency between interest rates and economic indi- 
cators, and serve as a basis for evaluating the role of non-equilibrium volatility in the bond 
market. The first answer we reach is that the equilibrium picture holds for an analysis of 
annually- averaged long-term interest rates. Looking more closely in a second analysis [8], we 
show that this picture breaks down for timeframes less than a year long during which non- 
equilibrium speculator behavior plays a key role, leading to large interest rate fluctuations. 
In particular, bond prices dropped from 57% to 21% of their face value, from July to Decem- 
ber 2011. A complex systems perspective suggests that near a transition large fluctuations 
between distinct states (in this case default and non-default) become possible. Under these 
conditions interest rate variations are both caused by, and cause, vulnerability to default. 
Such fluctuations can trigger substantial economic impacts and stronger interventions than 
would be justified by equilibrium markets. The extent of these fluctuations may be reduced 
by well chosen market regulations. 

II. OVERVIEW 

It is common practice for private corporations and governments to borrow in order to 
overcome shortfalls in operating expenses, deal with emergencies, or invest in economic 
growth. A primary mechanism for this borrowing is to sell interest bearing bonds. The 
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question of how much governments can borrow from bond markets is related to the lenders' 
concern about debt repayment and sustainability [9, 10]. The levels of government debt vary 
markedly from country to country, being the result of a mixture of policy decisions (if levels 
are low) and how much lenders will finance (if levels are high) . The interest rates are set by 
the willingness of those who loan funds, either by direct negotiation or by auction. Only if 
the issuer becomes unable to repay, i.e. a default happens, are the payments reduced. The 
interest on the primary bond market therefore is, in first approximation, a measure of how 
high the risk of default is. The more likely those who loan the funds think default is, the 
higher interest rates are set, reflecting a greater investment risk. Different bonds provide 
more or less interest depending on their estimated relative risk of default. 

Many researchers have addressed the question of why nations ever choose to pay off their 
debts, given sovereign immunity and no enforcing body to exact repayment [9-27]. Most 
assume that there is some cost associated with defaulting, which may include loss of reputa- 
tion in the international community [11-13], trade sanctions [16, 20], harsher future credit 
terms [14, 26, 28], or outright exclusion from the world credit market [9, 14]. Historically, 
military interventions may have served as an additional external deterrent [13, 17]. However, 
some question the importance of such external repercussions, claiming that they may not 
be as stringently enforced or as effective as traditionally thought [13, 15, 18]. By whatever 
mechanism, it is clear from measures of economic growth that sovereign default correlates 
with subsequent reduced economic performance of the defaulting country [13, 19-21, 29]. 

A complementary body of literature attempts to identify the warning signs that precede 
sovereign default [30-37]. Researchers have cited macroeconomic measures of insolvency 
and illiquidity as precursors to default; they consider GDP growth, debt-service payments, 
penalties for default, bond interest rates, interest volatility, consumption volatility, and 
a host of other factors [31, 33-35]. Some have validated their models on empirical data 
[31, 32, 34, 35, 38]. However, many of them include a large number of parameters, obscuring 
the relative importance of different factors. 

Here we develop a quantitative model of sovereign default that identifies the debt ratio [39] 
as the only relevant dynamic macroeconomic variable that European market participants 
have used over the last 10 years to estimate the likelihood of default of a country and, 
therefore, to set long-term interest rates. The model explicitly relates the annually-averaged 
long-term interest rates to the debt ratio with a simple 2-parameter fit, and is able to 
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accomplish two major goals. 

First, the equilibrium model provides a way of testing whether international capital mar- 
kets are consistently basing interest rates on economic conditions. We show that interest 
rates reflect the expectation that default becomes certain as the debt to GDP ratio ap- 
proaches a specific value, the default threshold. Although there is strong evidence for the 
existence of mechanisms able to drive prices away from equilibrium [40, 41], economic mod- 
els usually assume that markets are reflecting economic information rather than considering 
it a hypothesis to be tested. We test the assumptions of the equilibrium picture by showing 
that the annually-averaged long-term interest rates for Greek bonds over the last 10 years 
are given by a well-defined function of the distance of the sovereign debt ratio from a fixed 
country-specific default threshold. Interest rates follow this function toward default as the 
debt ratio increases (see Fig. 1). Since the interest rates can be quantified as an explicit 
function of macroeconomic information, the model shows that annually-averaged long-term 
interest rates are consistent with this assumption of equilibrium models. 

Second, the model can estimate the timing of a default event given its debt trajectory. 
Default becomes certain as the country debt approaches the debt threshold. Even when a 
country defaults a portion of the debt is typically honored, the amount of which is called 
the recovery rate, affecting the precise timing of the default. For Greece, we find that the 
model points to certain default with a 50% recovery rate in the first quarter of 2013 (Fig. 2, 
top). The "voluntary" write-down of half of the country's debt by private creditors, which in 
effect corresponds to a partial default with 50% recovery, occurred one year earlier in March 
2012 [42]. The timing of default may be attributed to non-equilibrium market fluctuations 
in the secondary bond market and not to economic or other news. Thus, the equilibrium 
trajectory of default is preempted by market bandwagon effects. Results for other countries 
are summarized in the bottom panel of Fig. 2. The projected timing of default ranges from 
March 2013 for Portugal to July 2016 for Italy. 

Our analysis shows, therefore, that the change of long-term interest rates follows funda- 
mental expectations only over the long run, while it is subject to non-equilibrium effects in 
the short run. We discuss this in a second paper, where we show that shorter time variation 
in bond prices have been much more rapid than would be expected from economic funda- 
mentals [8]. The results are consistent with the existence of price fluctuations due to trend 
following or market manipulation. These fluctuations are both caused by and cause vulner- 
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FIG. 1: Test of equilibrium market model of default - Greece's debt ratio Rt (blue dots) 
as a function of the parameter Xt, which reflects annually-averaged long-term interest rates. The 
equilibrium model derives a linear relationship between Rt and Xt (see Eq. 3). All points, from 
2001 to 2011, lie on a straight line, indicating that over this period of time markets have consistently 
determined a level of debt at which default occurs — the debt for Xt = (maroon star). 

ability to default. When there is a concern about default, large fluctuations are created by 
herd behavior. When interest rates fluctuate upwards they can cause default. Under these 
conditions the causal relationships are reversed and markets drive fundamental economic 
outcomes. These effects shorten the time to default, causing a need for more aggressive aus- 
terity programs and debt write-downs than might otherwise be required. During a period of 
economic recovery from a recession these differences may be critical as a recovering economy 
may be better able to avoid default if it is not subjected to higher national interest rates 
and austerity measures than would be justifled in equilibrium markets. 

In considering the broader implications of our analysis of country default, it is important 
to note that our model was developed and validated in the context of the European debt 
crisis, which is an atypical context due to the common euro currency that prevents individual 
country currency devaluation for members of the Eurozone. Where currency devaluation is 
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FIG. 2: Timing of default - Top Panel: Greece's debt ratio Rt normalized by our model param- 
eter Rc, the critical debt ratio, as a function of time (solid line). The projected debt trajectory 
(dashed line), intersects the normalized threshold in March 2013. Red mark indicates the actual 
partial default in March 2012. Bottom Panel: Similar analysis for all five EU countries considered 
in this work: Italy, Spain, Greece, Portugal and Ireland. See Appendix A for more details on the 
model, and Appendix D for more details on individual countries. 



possible, devaluation may serve to reduce some of the pressure of paying back internal 
debt, also by making exports cheaper and generating growth [43]. Thus, the European 
conditions may have distinctive properties and an analysis of European conditions may not 
be applicable elsewhere without modification. 

In particular, one of our central findings is that only one macroeconomic variable is 
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dynamically relevant, the debt ratio. The debt ratio has been recognized to be a common 
reference for lenders among measures that are expected to be important [10]. Our model 
allows for the possibility that other variables play a role in country default as long as they 
either do not change over time, and therefore contribute only to the static model parameters, 
or are not relevant in the specific context of the Eurozone countries we studied. Indeed, 
other variables must play a role in determining the value of the static country specific debt 
threshold. Our analysis does not reveal these variables. Still, the debt thresholds are in 
a limited range between 90% and 200% of GDP for the countries we studied — implying 
markets consider the GDP itself to be a first approximation to a reasonable debt threshold. 

This conclusion may, however, be modified for other countries outside the Eurozone. A 
variable that may not be relevant in the Eurozone because of the common currency, but 
could be relevant in other areas, is the distinction between external and internal debt [44]. 
Where one class of lender is substantially less likely to stop providing loans than others, 
an analysis may be best framed in terms of the distinct properties of the two levels of 
available debt. It is reasonable to expect that the internal debt plays a diminished role in 
the probability of default. If the ratio between the internal and external debt were constant 
over time, discounting the role of the internal debt in default would increase the value of 
the critical debt ratio, but would not otherwise affect our conclusions. 

Thus, we can expect our findings to be modified for countries where internal debt is a 
large part of total debt, especially for countries outside of the Eurozone. For example, Japan 
has a high debt ratio of 230% but much of the debt is internal and owned by the postal 
system [45]. Much of the US debt is owned by the Federal Reserve [46]. In both these 
cases, interest rates are low despite debt ratios that are in the range of the Eurozone debt 
thresholds. 

Finally, the applicability of the model is limited to countries in distress. We cannot know 
what the critical threshold of the debt ratio is when interest rates are low: the only way 
to test the model is when the probability of default is high enough to affect the annually- 
averaged long-term interest rates. 
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III. MODEL 



We construct a mathematical model of sovereign default risk (details are in Appendix 
A) by considering a country to have the choice, at each time step, of being in a state of 
default or nondefault. Since the ability of a country to pay its debt depends on the size of 
the debt relative to its economic output, we assume that the choice of default is determined 
by its debt ratio, i.e. the ratio between the overall debt and the GDP. We also assume the 
existence of a critical threshold Rc of the debt ratio, above which the country is likely to 
default. The value of the critical threshold depends on the cost of the penalties implicitly 
associated with defaulting; it differs from country to country and is the first parameter of 
our model. 

The expected probability of default Pt is indicated by market participants through the set- 
ting of interest rates: the more likely they think default is across all possible future scenarios 
over the period of repayment, the higher interest rates are set, reflecting a greater invest- 
ment risk. Since interest rates determine the debt burden, their increase also determines a 
higher probability of default. This positive feedback mechanism is shown as a blue dashed 
arrow in Fig. 3. Moreover, a sharp change in interest rates can alter investor perception of 
the default probability (orange arrow), thereby creating another feedback loop and pushing 
interest rates out of equilibrium. There are therefore two distinct feedback mechanisms, one 
fundamental and one behavioral. The model can test whether the behavioral feedback has 
an impact on the markets, or whether markets are self-consistently determining a default 
probability based on the fundamental feedback (blue arrows). 

If all lenders had the same expectation about the time of default given by a default 
threshold, the interest rate it would be the risk free value until the debt ratio Rt reaches 
the critical default value. At this point, creditors would stop lending money and interest 
rates would effectively diverge. The resulting step-like profile for the probability of default 
is in practice smooth, since their expectations are heterogeneous. Many reasons contribute 
to this heterogeneity, including the limited amount of information available to the economic 
actors in the bond market, different estimates of the critical debt ratio, imperfect knowledge 
of the debt trajectory, different evaluations of the political decision-making process behind 
a possible default, different bond maturation periods and payment schedules, and disparate 
influence of other economic indicators. Heterogeneous expectations have been previously 
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FIG. 3: Interest rate setting mechanism - In equilibrium, investors use macroeconomic indi- 
cators of national production and debt to estimate the probability of default (blue arrow marked 
1) using Eq. 1. The market then sets bond prices and interest rates according to that risk (blue 
arrow marked 2) using Eq. 6. Interest paid on the debt contributes to the subsequent debt (blue 
arrow marked 3), according to Eq. 22). The model can also determine whether sudden changes in 
interest rates driven by bandwagon effects increase the perception of default probability (orange 
arrow), pushing interest rates out of equilibrium. 

recognized to be important in describing the behavior of bond markets [47]. Nevertheless, 
the probability of default Pt is low below Rc and high above, so that the overall profile as a 
function of the debt ratio Rt can be represented by a smooth sigmoidal transition. 
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where ?7, the second parameter of our model, reflects the degree of smoothness, and p is the 
extrinsically set recovery rate, i.e. the probability that the country will pay a share of the 
debt p < 1. The derivation of this equation can be found in Appendix A. 

In summary, we consider the default process as one that follows a discrete jump at a par- 
ticular value of the debt ratio, but we treat it as smoothed by uncertainty and heterogeneity. 
The model has only two parameters: the critical debt ratio Rc-, which represents the average 
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FIG. 4: Default model validation: Greece - Plot of the increasing probability of default with 
increasing debt according to data for Greece (dots) and theory (dashed line, Eq. 1) using the 
relationship of interest rates to risk as specified by Eq. 2. The red star indicates the actual partial 
default in March 2012 [42], at a time when the probability of default was Pt = 0.33. The fit of 
theory to data yields the two model parameters, the critical debt ratio Rc = 2.00 it 0.07 and the 
heterogeneity parameter rj = 0.18 it 0.02. The data and theory shown are for 50% recovery. Vertical 
shading shows the range of debt ratios of certain default for a range of recovery rates between 20% 
to 80%. German interest rates are used as risk-free interest rates. Using quarterly data as opposed 
to annual data would introduce seasonality effects, but would not change significantly the results. 

value above which investors stop buying bonds, and the heterogeneity parameter 77, which 
incorporates market uncertainty. The model assumes that these parameters are well defined 
and consistently represented by the bond market over time. 

These two parameters can also be used to evaluate if the market interest rates of a 
country are consistently set according to economic fundamentals. The probability of default 
Pt increases as the interest rates it increase, and it goes to as they approach the risk-free 
interest rates r^, according to [50] 
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FIG. 5: Time of default: Greece - We obtain the time of default for Greece by projecting 
Greece's debt ratio (solid dots) as shown (dashed line). Default becomes certain when the debt 
to GDP ratio reaches the value Rc = 2.0 in March 2013. Default is preempted by the impact 
of the negotiations between Greece and the European Union including a 53.5% haircut, austerity 
and ESPF loans as shown in projections for 2012 and 2013 (green circles) according to estimates 
reported in Ref. [48]. The impact of deficit targets without the other interventions led to much 
higher projections that did not avoid default (blue circles) [49]. The horizontal shading corresponds 
to the range of debt ratios for defaults of 20-80% recovery (see Fig. 4). 

Equating Eq. 1 and 2 we have 

Rt = Rc- vXt (3) 

where Xt = ln(l-|-r() — ln(zi — r^) (see Appendix A for details). Therefore, if the experimental 
relationship between Rf and Xf is linear, international markets are consistently setting 
interest rates according to the debt ratio (see Fig. 1). 



IV. RESULTS 



We tested the model on Greece over the last decade. Results are shown in Fig. 4 assuming 
a recovery rate of 50%, the average value on defaulted sovereign bonds over the last 30 years 
[51]. (The model parameters do not depend on the recovery rate used.) The figure shows 
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the probability of default plotted against the debt ratio. As predicted by Eq. 1, the data 
lies on a sigmoidal curve, and the two model parameters can be extracted from this fit. We 
obtain a critical debt ratio Rc = 2.0, and a heterogeneity parameter 77 = 0.18. Therefore, the 
model estimates that when Greece's debt reaches a level about twice its GDP, the market 
projects the country default probability to be 1. 

Figure 5 shows how the time of certain default of a country can be estimated. We 
approximate Greece's debt trajectory in one of two ways, using a numerical polynomial fit 
or an economic debt projection (Appendix B). Either way, the debt ratio would reach the 
default threshold Rc in the first few months of 2013. Simulations with different intervention 
scenarios (Appendix D) show that the default could not have been prevented solely by 
meeting the austerity targets recently agreed upon [49]. Indeed, global fears of a Greek 
default led to negotiations and to a coordinated intervention that resulted in a haircut, 
which can be considered tantamount to the partial default projected by our model. The 
debt write-down occurred at a time when the probability of default was 33%, one year before 
the projected certain default. A separate analysis shows that the timing of default can be 
attributed to non-equilibrium market fluctuations [8]. Negotiated partial defaults are more 
likely than a full default because of the harm to both lenders and borrowers [52]. In this 
case, negotiations led to a 53% write-down of the debt by the private sector (in the form 
of a waiver of receivables) coupled with access to EFSF and IMF loans at preset interest 
rates and austerity measures. If the austerity measures are met, and its impacts on GDP 
are consistent with published expectations (Appendix C), this will allow the debt ratio to 
remain below the danger zone. 

We can use our results to evaluate whether market interest rates are consistent with 
equilibrium assumptions over time. During the near-default period, bond rates change 
rapidly and are subject to higher uncertainty. In this scenario, the impact of large traders 
and trend following may become relevant. The fit of the default dynamics for Greece shown 
in Figures 4 and 7 is a test of market self-consistency: if the market did not set interest rates 
consistently with the valuation of the critical debt ratio in prior years, these relationships 
would not hold over time. For example, if the interest rates overestimated the probability 
of default of a country in a given year relative to the historical precedent, the corresponding 
point would lie above the fit in Fig. 4. The self-consistency condition can also be represented 
as a linear relationship between debt ratio and the distance of market interest rates from 
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their default value (see Fig. 7). The good fits imply that the reference 10 year bond interest 
rates averaged over yearly time frames are consistent with market equilibrium. Shorter time 
frame behaviors deviate from these conclusions and are treated in a separate paper [8]. 

We also tested the model for the other European countries involved in the debt crisis: 
Portugal, Ireland, Spain and Italy (see Fig. 9-12 in Appendix D). As in the case of Greece, 
we found that markets self-consistently determined a definite value for the critical debt ratio 
Rc for all countries. The time range over which a fit with a consistent value of Rc is possible 
depends on the country: the model fits the last 11 years for Greece and Italy, 9 years for 
Portugal, and 5 years for Spain and Ireland. 

As for the case of Greece, the equilibrium model portrays a situation which could be less 
severe than the one described by secondary markets and media [53-56] . If the budget targets 
are met, Portugal could avoid a partial default at the beginning of 2013. Ireland and Spain 
have a larger buffer before a default, projected for mid-2014. Italy is still farther away from 
a possible default, projected by the model for mid-2016, if the current trend is maintained. 
In all cases an immediate large scale policy response to change the debt trajectory does not 
seem to be necessary. 

The values of the critical debt ratio, Rc, do not vary over orders of magnitude between 
countries. The values cluster around unity, which could imply that lenders believe a country 
is able to repay its debt when the debt is about the size of its GDP, but not significantly 
larger. The capacity of a country to borrow depends on a variety of factors, and reflects 
creditors' concerns about solvency (long-run prospects), liquidity (short-run prospects) and 
repudiation risks [10]. The reasons why countries usually borrow lower-than-sustainable 
amounts, i.e. why Rt is kept well below Rc, have been discussed [10]. Even though we find 
there is only one scale for this parameter, Rc varies among countries. The value of Rc may 
depend on a variety of macroeconomic phenomena. Our analysis does not identify these 
dependencies. Sustainable growth rate, the coupling to other economies willing to help, 
considerations about the country's hidden economy and its possibility to increase taxes are 
also factors that might contribute to the country-specific value of Rc- 

Finally, it is worth noting that the framework we presented here is a long-run model 
of default for European countries, but short timescales require a different picture. This is 
demonstrated by the deviation of 10-year Greek bond prices from their equilibrium trajectory 
after the summer of 2011, when they dropped from 57% to 21% of their face value in five 
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months. One might think that this increase in interest rates would be justified by the default 
projections, even if they were about 1.5 years away according to our equilibrium model. But 
the corresponding probability of default at that time was around 0.3 (see Fig. 4), not large 
enough to explain such a high level of interest rates, nor to justify an imminent write-down. 
The equilibrium economics perspective is uncompromising: bond prices should follow the 
curve regardless of fluctuations — the trajectory has been deflned. In contrast, a complex 
systems perspective predicts that once the system is near a critical point, multiple stable 
states (in this case default and non-default) become possible. In this regime, the market 
is not self-averaging, its flnite size gives rise to fluctuations at all scales, and sensitivity to 
trend following and market manipulation increases. Both equilibrium and non-equilibrium 
mechanisms are therefore included in our framework and are not inconsistent, as they are 
acting at different timescales. 



V. CONCLUSIONS 



We developed a quantitative model of sovereign default for European countries that 
identifies the debt ratio as the fundamental parameter that investors use to set the persistent 
value of interest rates captured by annually-averaged values. Our analysis shows that the 
equilibrium model assumption that interest rates are set by economic conditions is satisfied, 
though we do not show that interest rates actually refiect potential losses due to default 
risk. Our model is able to address several questions. First, it gives an estimate of the debt 
sustainability of a country, i.e. how much creditors are willing to lend. In the case of Greece 
the maximum debt is twice the GDP, R^. ~ 2.0. Second, given the debt trajectory of the 
country, the model can identify the time the debt ratio goes beyond the threshold, triggering 
a default. Since default is itself determined by the time of diverging market interest rates, 
default is self-consistently determined by the market evaluation of the default probability. 
Third, the model can establish whether bond markets satisfy the equilibrium assumption 
that interest rates the market sets have a consistent relationship with fundamental indicators 
over time. We find this to be the case for Greece over the last 10 years. 

The model identifies the time of certain default as occurring early in 2013. While this 
appears to be in reasonable agreement with the debt write-down, the timing is different by 
approximately a year and this is critical to a proper interpretation of events. The existence of 
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non-equilibrium bandwagon effects that preempt the equihbrium default means that higher 
national interest rates and austerity measures than would be justified in equilibrium accel- 
erated and exacerbated the default process. Without these effects, more modest measures 
and less social disruption could have been sufficient to avert default. 

As it stands today, drastic measures taken by the government, international organiza- 
tions and the private sector may avoid a complete default. These interventions have both 
an economic and social cost. By regulating markets to limit the non-equilibrium market 
fiuctuations, the cost of interventions and their severity could be reduced for Greece and 
other European countries. 

VI. ACKNOWLEDGEMENTS 

We thank Karla Z. Bertrand for help with the literature review, Jeffrey Fuhrer, Richard 
Cooper, Yavni Bar- Yam and Dominic Albino for helpful comments. 



16 



Appendix A 



Model Details 



We consider a country that decides to borrow from international capital markets at 
interest rate i. At every time period t, the country is in one of two states: default or 
nondefault. Since the ability of the country to pay its debt D depends on its economic 
output Y, the state of the country is determined by its debt ratio, R = D/Y. When R 
is above a critical threshold Rc, the value of which depends on the cost of the penalties 
associated with defaulting, the country is in the default state, and vice versa. 

If P is the probability of default with recovery rate p, i.e. the probability that the country 
will pay a share of the debt p < 1, the expected return from a one-period loan would be 
{1 + — P + Pp) [50]. By definition, (1 + i) is the contractual repayment at the end of the 
period per unit of loan and (1 — P) is the probability of full repayment. An investor, however, 
has the option of investing in risk-free bond markets and receive an expected return of (1 + r) 
at the end of the period, where r is the risk-free interest rate. Thus, since in competitive 
equilibrium investors should be indifferent between those two alternatives, equating the two 
expressions for expected return and solving for P we have 



This equation establishes a relationship between the expected probability of default and 
interest rates. It is consistent with the intuition that the probability of default P increases 
as the interest rates i increase, and it goes to as they approach the risk-free interest rates. 
Equation 4 can also be written as a function of bond prices, B. For a one-period loan, 



where Bq is the initial value of the bond (face value) and io the initial interest rate. Then, 
substituting in Eq. 4 we get 




(4) 



B 



l + i 



(5) 




(6) 



which relates the probability of default with the current bond price. 
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There are three reasons why i is in general greater than r: transaction costs, risk aversion 
and default risk [50]. For simplicity, let's assume that transaction costs are small enough 
to be neglected and that lenders are risk neutral. If lenders had perfect knowledge of the 
time of default of the country, and the duration of the one-period loan is short enough to be 
neglected, they could keep i = r until the debt reaches the critical default value Rc. At the 
time of default, when the debt reached the critical default value, the interest rates would 
become infinite. Therefore, at each time t this idealized bistable system can be summarized 
step function of the default probability: 

Pt = 0, it = rt Rt < Rc 

(7) 

Pt = l, it = oo Rt> Rc 

We plot the default probability for the case Rc = I 'm. Fig. 6 (black curve). This step func- 
tion is in practice smooth (blue curves), since the investment strategies are heterogeneous for 
a number of reasons: limited amount of information available to market participants, differ- 
ent estimates of the critical debt ratio, imperfect knowledge of the debt trajectory, different 
evaluations of the political decision-making process behind a possible default, different bond 
maturation periods and payment schedules, disparate influence of other economic indicators. 
Here we show how this smooth behavior naturally arises from a supply and demand model 
of capital markets that includes heterogeneous expectations. The result follows from the 
reasonable assumption that all dependencies are smooth other than the singular behavior of 
the default itself. Smoothness in the vicinity of the default justifies expanding all functional 
relationships to first order. The following derivation is illustrative and captures the eventual 
behavior. 

If we assume a linear relationship between bond prices B and quantity of debt supplied 
by the country {Qs) or demanded by investors {Qd) we have 

Qs = as + /3sB (8) 
Qd = ad- (5dB (9) 

and, since at equilibrium Qg = Qd, 
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(10) 



assuming that supply and demand shocks would influence the intercepts — and not the 
slopes — of the linear dependences. The demand intercept ad depends on the overall liq- 
uidity available for financing the country's debt, and therefore on the number of market 
participants. Rational investors are willing to lend money to the country as long as its debt 
ratio is below their estimate of Rc- Due to all of the factors contributing to heterogeneity, 
we expect investor estimates of the critical value of default to have a distribution, whose 
form can be approximated by a normal distribution with average Rc and standard deviation 



f{Rt) 



N 



exp 



1 ( Rt - Rc 

2 V o 



(11) 



where iV is the total number of investors. Since agents that estimate the default at a 
particular value of Rt expect it also for any greater value of the debt ratio, is proportional 
to the cumulative distribution function 



N 
~2 



erf 



Rt — Rc 



a 



(12) 



where erf () is the error function. This expression can be well approximated with a sigmoid 
function [57], according to 



N 



(13) 



1 + (,{Rt-Rc)/'n 

where rf = crA/27r/4 is the parameter characterizing the heterogeneity. The bond supply 
intercept, as, depends on the amount of money a country wants to borrow, and thus on 
its debt ratio. We expand this dependence around the critical value Rc, keeping only the 
constant term. The effect of the absolute level of the debt ratio (-Rt) on the supply (and more 
generally on bond prices) can be neglected compared to the effect of the sigmoidal function 
dependence {Rt — Rc), i.e the distance of the debt ratio from its critical value. Substituting 
Eq. 13 into Eq. 10 we have 



Bt = k' 



1 + (.^Rt-R.)h 



+ k" 



(14) 
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FIG. 6: Default probability - Dynamics of the probability of sovereign default as a function 
of the country's debt ratio, assuming no recovery. In an idealized scenario in which lenders have 
perfect information, a step function of the default probability would occur (black line, Eq. 7). 
Imperfect knowledge and heterogeneous strategies results in smoothed sigmoidal function (blue 
lines, Eq. 16). The extent of smoothing depends on the width of the distribution of expectations 
given by a. 

where k' and k" are constants that are functions of a^, /3s, i^d and iV. Their value can be 
inferred by imposing constraints on two hmiting cases of Eq. 5: 

• it — oo, -Bf — )■ when i?t — )■ oo 

• it^ n, Bt Bo{l + io)/{l + rt) when Rt ^ 



Using these conditions we obtain 



Bt = Bq 



l + io\ l + e--f^=/^ 



Bo 



1 



.l + nJ 1 + e^^^-^-)h V 1 + n / 1 + e'^^^-^-)l'^ 

since e~^''f^ <^ 1. Inserting this expression in Eq. 6, 



(15) 



Pi 



1 



1 + (,{Rt-Rc)h 



1 



l-p 



(16) 
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FIG. 7: Test of equilibrium market model of default - Greece's debt ratio [3] as a function 
of the default distance parameter Xt (Eq. 19) from 2001 to 2011 inclusive. Xt reflects long-term 
interest rates, i.e. secondary market annually-averaged yields of 10-year bonds [59]. Data (solid 
line) can be fitted with Eq. 19 (dashed line) to obtain the two model parameters, Rc = 1.90 it 0.05 
and r] = 0.15 it 0.01. The small difference between these values and the ones obtained from the 
fitting in Fig. 4 is due to the logarithmic scale. Data shown for the debt ratio are end-of-year 
estimates, while long-term interest rates are averaged over the course of the year. Using annually- 
averaged debt ratios or quarterly values for both debt ratio and interest rates does not significantly 
affect the results. Fitted parameters in all cases are the same within statistical uncertainty. 

This equation relates the probability of default to the two parameters of our model, Rc and 
1]. The family of curves Pt vs fZf at p = for different values of a are shown in Fig. 6. The 
asymptotic behavior of the probability of default approaching 1 when the debt ratio goes to 
infinity is justified by theoretical [9], empirical [58] and numerical data [38]. 
Equating Eq. 4 with Eq. 16 gives 

R, = R^-^ln(^±^) (17) 
V«t -rtJ 

When the debt ratio reaches the critical threshold, Rt = Rc, the corresponding critical value 
for interest rates is 

= 1 + 2ri ~ 1 (18) 
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since risk-free interest rates are negligible compared to the interest rates of a country close 
to default. We can rewrite Eq. 17 as, 

Rt = Rc- TlXt (19) 

where Xt is a measure of the distance between the current value of interest rates and the 
value at default. This linear relationship between Rt and Xt is verified empirically in Fig. 
7 for Greece, and in Appendix D for Portugal, Ireland, Spain and Italy. While the slope of 
the linear fit gives the heterogeneity parameter, its intercept is the critical debt ratio of the 
country. 

The debt threshold Rc does not depend on the assumed amount of debt recovery, but 
the time of expected default does. Eq. 16 specifies the probability of default for different 
debt thresholds and implicitly the debt ratio at which the default probability reaches one, 
Rd = Rc + ?7ln(l/p — 1). For the generally observed case of 50% recovery, the critical debt 
ratio coincides with the debt threshold, Rd = Rc. 
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Appendix B 
Debt Feedback Loop 



In this Appendix, we provide a quantitative description of the feedback loop in Fig. 3 
(blue arrows). We are going to use this analysis to identify the expected change of the 
debt ratio over time, given policy options that provide below market interest rates such as 
austerity and bailouts. 

The country's debt trajectory is 

Dt = St + Dt-,{tt + l) (20) 

where D is the gross debt, i is the market interest rates and S is the primary deficit, i.e. 
government spending minus tax revenue. If we divide the previous equation by the economic 
output Yf we have 

Rt = st + ^itt + l) (21) 

where St = St/Yf is the budget ratio. If we define the GDP growth as yt = (Yt — Yi-i)/Vt-i 
and note that Dt-i/Yt = Rt-i/{l + yt), from the previous equation we have the accumulation 
equation, 

Rt = st + Rt-i ( + l) . (22) 

This equation completes the feedback loop represented in Fig. 3: fundamental economic 
indicators determine the expected default probability (arrow 1 and Eq. 16), which influences 
market interest rates (arrow 2 and Eq. 6), which in turn influences the new value of economic 
indicators (arrow 3 and Eq. 22). This feedback loop reflects the role of interest payments 
themselves on the increase in country debt. 

Finally, we show how we can solve Eq. 22 and fit the debt trajectory analytically. The 
following equations, which are not used anywhere else in the paper or in the default model, 
demonstrate that in order to avoid default in the not-too-long term, the GDP has to grow 
faster than borrowing. As we saw in Eq. 22, interest rates can affect the debt ratio, and 
the debt ratio grows very rapidly if it > yt, i.e. if the interest rates are larger than the rate 
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FIG. 8: Extrapolation - Greece debt ratio trajectory projection using Eq. 25 and assuming 
time-independent quantities, with quarterly values: i = 0.021, s = 0.008 and y = —0.007. 

of economic growth. To obtain a closed-form solution that shows the dependence of the 
debt ratio on the GDP change and interest rates, we assume that it, yt and St are not time 
dependent — this assumption is only used for the analytic fit in this Appendix, and not for 
our model or elsewhere in this paper — and solve Eq. 22 to give 



Rt — Ro 



1+zV 1+y 
+ s- 



1 



(23) 



\l + yj i-y 

where Rq is the value of the debt ratio at time 0. The difference equation 22 can be 
transformed into a differential equation if the usual assumption of smoothness applies. 



dR{t) , y 
dt 



l+y 



Rit) 



(24) 



This equation can be solved to give 



R[t) = Roexp I :;— — t I + s 



l + y 



y 



I — y ' 
exp 1 



(25) 
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Equation 23 coincides with Eq. 25 for small values oii — y, and either demonstrates how the 
debt ratio grows exponentially iii > y. Since our model shows that there is a threshold debt 
ratio set by markets, an exponential growth is unsustainable. The fit of the time dependence 
of the Greek debt ratio with Eq. 25 is shown in Fig. 8. We emphasize that Equation 22, 
which includes the time dependence of it-, yt and St, is more general than Equation 25. The 
latter is only used in this Appendix and not elsewhere in this paper. 



25 



Appendix C 
Austerity Measures 

In order to evaluate the effect of austerity measures on tlie debt ratio of a country, we 
use Eq. 22 to calculate Rt given estimates by economic analysts for the GDP change yt and 
austerity targets for the budget ratio st. We assume the interest rates for the next two years 
are similar to current rates, it ~ it-i-, which is justified by the small absolute difference 
over time of the short-term interest rates, and the very few years over which the projection 
is being calculated. Our purpose is solely to evaluate whether the magnitude of austerity 
measures that have been implemented may be able to avert a default, using available eco- 
nomic projections according to the debt threshold we obtained. We report the estimates for 
yt and the targets for st in Table I. The impact of the austerity programs on debt ratios are 
incorporated into the figures in Appendix D, showing their implications for averting defaults. 



Country 


"S2012 


S2013 


1/2012 


Z/2013 


bailout 


Greece 


7.3 [49] 


4.7 [49] 


-6.4 [48] 


-1.9 [48] 


yes 


Portugal 


4.5 [60] 


3.0 [60] 


-3.4 [61] 


0.0 [61] 


yes 


Ireland 


8.6 [62] 


3.0 [63] 


0.7 [64] 


2.2 [65] 


yes 


Spain 


5.3 [66] 


3.0 [66] 


-3.0 [67] 


-1.9 [67] 


no 


Italy 


1.7 [68] 


0.5 [69] 


-1.2 [70] 


0.5 [69] 


no 



TABLE I: Values used for debt ratio projections - Austerity targets for the budget ratio st 
and projections of the GDP change yt for the five countries we considered. 
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Appendix D 
Other Countries 



A. Portugal 
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FIG. 9: Model fits and projections for Portugal - Top left •panel: Portugal's debt ratio [3] as 
a function of the default distance parameter Xt (Eq. 19) from 2001 to 2011 inclusive. reflects 
long-term interest rates, i.e. secondary market yields of 10-year bonds [59]. Data (solid line) 
can be fitted with Eq. 19 (dashed line) from 2007 to 2011 to obtain the two model parameters, 
i?c = 1-35 lb 0.05 and r/ = 0.11 it 0.01. Top right panel: Probability of default as a function 
of Portugal's debt ratio. Pt is calculated from Eq. 4, setting the recovery rate to p = 0.5 and 
considering German interest rates as the risk-free interest rates. Data (dots) can be fitted with 
Eq. 16 (dashed line). This allows us to extract the two model parameters, the critical debt ratio 
Rc = 1.34 lb 0.07 and the heterogeneity parameter ry = 0.10 it 0.02. Bottom panel: Portugal's 
debt ratio (solid dots) as a function of time. The debt trajectory, projected with a third-order 
polynomial regression (dashed line), intersects the critical debt ratio {Rc = 1.34) in Mar 2013. 
Hollow circles are projected debt ratios for the end of 2012 and 2013 if austerity measures are met 
(calculated from [60, 61]). 
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B. Ireland 




FIG. 10: Model fits and projections for Ireland - Top left panel: Ireland's debt ratio [3] as 
a function of the default distance parameter Xt (Eq. 19) from 2001 to 2011 inclusive. Xt reflects 
long-term interest rates, i.e. secondary market yields of 10-year bonds [59]. Data (solid line) 
can be fitted with Eq. 19 (dashed line) from 2007 to 2011 to obtain the two model parameters, 
Rc = 1.57 lb 0.14 and 77 = 0.20 it 0.03. Top right panel: Probability of default as a function of 
Ireland's debt ratio. Pt is calculated from Eq. 4, setting the recovery rate to p = 0.5 and considering 
German interest rates as the risk-free interest rates. Data (dots) can be fitted with Eq. 16 (dashed 
line). This allows us to extract the two model parameters, the critical debt ratio Rc = 1.63 it 0.10 
and the heterogeneity parameter r] = 0.21 it 0.03. Bottom panel: Ireland's debt ratio (solid dots) as 
a function of time. The debt trajectory, projected with a linear regression (dashed line), intersects 
the critical debt ratio {Rc = 1.63) in Apr 2014. Hollow circles are projected debt ratios for the end 
of 2012 and 2013 if austerity measures are met (calculated from [62-65]). 
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C. Spain 




FIG. 11: Model fits and projections for Spain - Top left panel: Spain's debt ratio [3] as a 
function of the default distance parameter Xt (Eq. 19) from 2001 to 2011 inclusive. Xt reflects 
long-term interest rates, i.e. secondary market yields of 10-year bonds [59]. Data (solid line) 
can be fitted with Eq. 19 (dashed line) from 2007 to 2011 to obtain the two model parameters, 
Rc = 1.02 lb 0.08 and rj = 0.10 ± 0.01. Top right panel: Probability of default as a function of 
Spain's debt ratio. Pt is calculated from Eq. 4, setting the recovery rate to p = 0.5 and considering 
German interest rates as the risk-free interest rates. Data (dots) can be fitted with Eq. 16 (dashed 
line). This allows us to extract the two model parameters, the critical debt ratio Rc = 0.93 ± 0.05 
and the heterogeneity parameter rj = 0.07 it 0.01. Bottom panel: Spain's debt ratio (solid dots) as 
a function of time. The debt trajectory, projected with a linear regression (dashed line), intersects 
the critical debt ratio (Rc = 0.93) in May 2014. Hollow circles are projected debt ratios for the 
end of 2012 and 2013 if austerity measures are met (calculated from [66, 67]). 
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D. Italy 




FIG. 12: Model fits and projections for Italy - Top left panel: Italy's debt ratio [3] as a 
function of the default distance parameter Xt (Eq. 19) from 2001 to 2011 inclusive. Xt reflects 
long-term interest rates, i.e. secondary market yields of 10-year bonds [59]. Data (solid line) 
can be fitted with Eq. 19 (dashed line) from 2007 to 2011 to obtain the two model parameters, 
Rc = 1.51 lb 0.1 and ij = 0.08 it 0.02. Top right panel: Probability of default as a function of 
Italy's debt ratio. Pt is calculated from Eq. 4, setting the recovery rate to p = 0.5 and considering 
German interest rates as the risk-free interest rates. Data (dots) can be fitted with Eq. 16 (dashed 
line). This allows us to extract the two model parameters, the critical debt ratio Rc = 1.44 ±0.04 
and the heterogeneity parameter rj = 0.06 it 0.01. Bottom panel: Italy's debt ratio (solid dots) as 
a function of time. The debt trajectory, projected with a linear regression (dashed line), intersects 
the critical debt ratio {Rc = 1.44) in Jul 2016. Hollow circles are projected debt ratios for the end 
of 2012 and 2013 if austerity measures are met (calculated from [68-70]). 
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